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ABSTRACT 

We show how certain properties of the Anderson model on a tree are related to the 



■ solutions of a non-linear integral equation. Whether the wave function is extended or 
^ ! localized, for example, corresponds to whether or not the equation has a complex solution. 
^ ' We show how the equation can be solved in a weak disorder expansion. We find that. 



for small disorder strength A, there is an energy E(.{X) above which the density of states 
and the conducting properties vanish to all orders in perturbation theory. We compute 
S ' perturbatively the position of the line Ec{X) which begins, in the limit of zero disorder, 
at the band edge of the pure system. Inside the band of the pure system the density of 
states and conducting properties can be computed perturbatively. This expansion breaks 
down near Ec{X) because of small denominators. We show how it can be resummed by 
choosing the appropriate scaling of the energy. For energies greater than Ec{X) we show 
that non-perturbative effects contribute to the density of states but have been unable tell 
whether they also contribute to the conducting properties. 
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1. Introduction 

Anderson localisation [1] , or the study of transport properties of a quantum particle in a 
random potential, is one of the most important problems in the theory of disordered systems 
[2,3]. In one and two dimensions an arbitrarily small random potential suffices to localise 
all energy eigenstates. In three and higher dimensions both localised and extended states 
can exist: strong disorder or energies far from the band center give rise to localised states 
whereas weak disorder and energies close to the band center produce extended states. 
Extended and localised states are separated by a line in the energy-strength of disorder 
plane, the mobility edge. The location of the mobility edge is a question of fundamental 
interest [4]. 

As usual in statistical mechanics, the simplest cases one can consider are mean field 
models. The most extensively studied mean field model of localisation is the Anderson 
model on a tree [5-12]. Various approaches have been developped, based in particular on 
self energy calculations [4,5] or on supcrsymmctry [8,9,13], which reduce the problem to 
a non-linear integral equation [5-9]. This integral equation, however, is complicated and 
the position of the mobility edge cannot be determined without recourse to some kind of 
approximation. Several works tried to overcome this difficulty by considering simplified 
versions of the model on a tree [14,15]. 

In the present paper, we reconsider the Anderson problem on a tree. We first give a 
derivation of the integral equation to be solved which, although completely equivalent to, 
is, we think, more intuitive than previous derivations. The system insulates or conducts 
depending on whether the integral equation possesses real or complex solutions. We try 
to solve this equation in the limit of weak disorder using a method [16] which generalises 
previous weak disorder calculations in one dimension [17]. 

One interesting outcome of this approach is the existence of a line Ec{\) in the E, A 
plane {E is the energy and A measures the strength of disorder) beyond which the integrated 
density of states and the conducting properties vanish to all orders in perturbation theory. 
This line tends to the band edge of the pure system, E = 2\/K (where K + 1 is the 
coordination number of the tree) in the hmit of zero disorder. We can show that for 
energies greater than Ec{X), non-per turbative contributions to the density of states make 
it non-zero. We have not, however, been able to determine whether non-perturbative effects 
also contribute to the conducting properties. The question is of particular interest because 
Abou-Chacra and Thouless [5] predict that the mobility edge tends to E — K + 1 rather 
than to the band edge, E — 2y/K in the limit of zero disorder. 

The paper is organised as follows: In section 2, we derive the non-linear integral 
equation satisfied by the distribution P{R) of a Riccatti variable R, defined to be the ratio 
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of the wave function at adjacent sites on the lattice. We show how the solution of this 
non-linear integral equation gives the integrated density of states and why the existence of 
a complex solution is related to the existence of extended states. In section 3, we discuss 
the pure system, i.e. the problem in absence of disorder. In section 4, we show how a weak 
disorder expansion can be performed for energies inside the band of the pure system. We 
find that in the presence of weak disorder the system conducts. In section 5, we extend 
the weak disorder expansion to the neighborhood of the band edge. We obtain within this 
perturbative approach an expression for the mobility edge Ec{X) in powers of the strength 
A of the disorder. For E > Ec{\), the integrated density of states vansihes to all orders in 
A, although it is known that for distributions of the potential with unbou nded support, 
it never vanishes [18-20]. In section 6, we discuss the origin of non-perturbative effects for 
energies outside the band of the pure system. Lastly, in section 7, we describe a numerical 
method to obtain the mobility edge, and we compare the results of this approach with the 
prediction of section 5 and with an exactly soluble case where the Ricatti variables are 
independent. 

2. Formulation of the Problem 

We consider a tight binding model on a Cayley tree of sites (see figure 1) with a 
random potential Vi at each site i of the lattice. The potentials Vi are independent random 
variables governed by a probability distribution p{V) which we choose to have zero mean 
{{Vi) = 0). The Schrodinger equation reads 

K+l 

J2 i^j = Eiji-XVi^i. (1) 

3=1 

Here ipi is the value of the wavefunction at site i, A is a parameter that controls the strength 
of the random potential, E is the energy of the particle, and the sum is over the K + l 
neighbors of the site i. It is useful to rewrite (|l|) as a recursion relation [16]. Call the 
central site of the tree io and define a Ricatti variable Ri on a site ihj Ri = tpj/tpi, where j 
is the neighbor of site i closer to io on the tree. Dividing (|^) through by ipi and regrouping 
terms gives (figure 1) 

R^ = E-XV-j2^ (2) 

j=i 

for all sites except io- This recursion allows one to calculate the Ri associated to all sites 
of the tree except for the site io (where it is not defined) and except for the sites adjacent 
to the boundary, where the Ri depend on the boundary conditions and, as we will see 
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later, should be chosen differently depending on the properties we want to study (density 
of states or conducting properties). 

On account of the random potential in (|^), the Ri are random variables governed 
by a probability distribution. The recursion (H) completely determines the probability 
distribution Pm{Ri) of an Ri located m steps from the boundary of the tree once the 
probability distribution Pi(-Rj) of the Ri on sites adjacent to the boundary have been 
specified. In what follows we will always choose boundary conditions in such a way that 
the Ri on sites adjacent to the boundary are identically distributed with distribution Pi{R). 
The symmetry of the tree then ensures that all the Ri an equal number m of steps from 
the boundary are also identically distributed with a probability distribution Pm{Ri) that 
depends on the number m one has to iterate (|^). The recursion (Q) induces the following 
recursion on the distributions Pm{R) 

Pm+iiR)= I Y[PmiRj)dRj J p{V)dV 5(^R-E + XV + J2^^ . (3) 

We shall assume in what follows that, for all the boundary conditions we consider {Pi{R) 
concentrated on the real axis when we calculate the density of states or Pi{R) concentrated 
on the lower half complex plane with all the Ri having a negative imaginary part when 
we study the conduction properties), the recursion @ converges to a limiting distribution 
P{R) which satisfies 

P{R) = I n PiRj)dR, J p{V)dV ^{R-E + XV + jz^^ . (4) 

Up to a change of variables, this integral equation is equivalent to the integral equations 
obtained in [5-9]. A similar equation also exists for diluted lattices [21-23]. The particular 
limiting distribution to which converges might depend on the initial Pi{R). We will see 
below that the localised and the extended regimes correspond to one of the two following 
situations: 

In the localised region: 

There is only one fixed distribution, Pi.eai(-R), which solves (|). This distribution 
is concentrated on the real axis and is stable, i.e. it is the limit of the sequence 
Pm{R) obtained through the recursion for any initial distribution Pi{R) (all initial 
distributions concentrated on the real axis as well as those concentrated initially in 
the complex plane converge to this distribution Pi.eai(-R), so that even if the R have 
initially some imaginary part, they become real under the iteration of (|^)). 
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In the extended region: 

There exist two different fixed distributions, Pj-eaiiR) and -Pcompiex(-R), which solve 
@. The real distribution Preai(-R), concentrated on the real axis, is the limit of 
the sequence Pm{R) when the initial distribution Pi{R) is concentrated on the real 
axis. This distribution Preai(-R) is however unstable against imaginary perturbations: 
a small imaginary component in the Ri on the boundary will not vanish under 
iteration of @. Instead, if, as for the scattering situation described below, the initial 
distribution Pi (R) is concentrated in the lower half complex plane (so that the initial 
R have all a negative imaginary part), Pm{R) converges to a different distribution 
-fcompiex(-R) concentrated in the lower half plane. (Note that there exists also a third 
distribution in the upper half plane symmetric to Pcompiex(-R) but we won't consider 
it because with the boundary conditions we use all the Ri are always either real or 
complex with negative imaginary parts.) 

The two fixed distributions Preai(-R) and Pcompiex(-R) are both solutions of the fixed point 
equation (||) and a great deal of what follows is devoted to the study of these fixed 
distributions. 

For choices of E, A and p{V) such that (^ has both a real and a complex fixed 
distribution (the extended phase), Preai(-R) and Pcompiex(-R) are not independent. In the 
appendix we show that the real fixed distribution Preai(-R) solution of (^) is given in terms 

of Pcomplex(^) by 

]^ poo roo ^ 
-Preal(-R) = " / dr ds— -5- ^-Pcomplcx (?^ " is) 

TT J-00 JO (P — r)^ + s'' 

— I r roo roo I ^ 

= — lm< dr ds— —PcompieAr-is)\. (5) 

71 U-oo Jo R — r + is J 



We now address the question of the choice of the initial distribution Pi(P) and discuss 
how the fixed distributions Preai(-R) and Pcompiex(-R) are related to the density of states and 
to the conducting properties. 

The density of states: 

Let us first discuss how the density of states can be obtained for the tree geometry [16] 
(see figure 1). We want to calculate the eigenvalues with the boundary condition that the 
wavefunction vanishes on the boundary of the tree. With this boundary condition, the 
Schrodinger equation for a site i adjacent to the boundary reads 

iPj = E^, - XV.^jji (6) 
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where j is the only neighbor of site i on the tree. Dividing through by ipi then gives 

R, = E-XV,. (7) 

where Ri = ipj/ipt- The initial R are real, so they remain real under iteration of (@) and 
the invariant measure P{R) is concentrated on the real axis. 

As discussed in [16], the equations (H) and solve the Schrodinger equation 
everywhere on the tree except on the central site iq. In terms of the Ri, the Schrodinger 
equation for the central site iq reads 

K+l I 

E-XV,,-J2 — =0 (8) 

j=i 

where the sum runs over the K + l neighbors j of site i^. All the Ri are functions of E and 
the values Ea of the energy which satisfy (H) are the eigenenergies. 

Expression contains all the information on the density of states but is not easy to 
use. If, however, one multiplies (||) by the product of all the Ri in the lattice [16], it becomes 
a polynomial in E of degree where A^ is the number of lattice sites and the coefficient of 
highest degree is 1 (for a tree of depth n, the number of sites A^ = [K"'-'rK^~^ — 2) / (K — l)). 
Therefore, one can write 

N / K+l 1 \ 

mE-Ea)=\E~XV,,-j:-^] I[R^ (9) 

Q = l \ j = l / i^io 

Both sides of this equation are polynomials in E with real coefficients and real roots. To 
extract the density of states, one can take the logarithm of this equality, for any complex 
value of the energy E, with the convention that the branch cut runs along the real axis from 
— oo to the largest eigenvalue Ea- When the energy approaches the real axis at a certain 
value E from above, the imaginary part of the left hand side is just vr times the number of 
eigenenergies larger than E whereas the imaginary part of the right hand side is equal to 
TT times the number of negative Ri (+1 when the term E — XVi^ — J2j is negative). This 
is because all the Ri as well as the term E — XVi^ — J2j have positive imaginary parts 
when the energy E is in the upper complex plane (see (0) and(^). Therefore the number 
fln{E) of eigenvalues greater than E (the integrated density of states) is given by 

K+l 1 

n^{E) = Qixv,, + y^—-e) + y: Q{-R^) (10) 

j=l i i^io 

Since the number of negative Ri is equal to the number of nodes of the wave function, we 
see that the equality ( pUj ) between the integrated density of states and the number of nodes 
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of the wave function, well known in one dimension, remains valid for tree structures as was 
already discovered by Dhar and Ramaswamy [24] in a similar calculation of the eigenmodes 
of Eden trees. 

For a tree of depth n (with N = + K"- — 2)/{K — 1) sites), the average over 

disorder of Qn{E) is given by 




+ + 1) K^-"^ f PUR)dR . (11) 

m=l •'-'^ 

This expression can be simplified by using the recursion (^ 

Pn{R) dR / Pn+i{R') dR! Q{ -B!) + {K+l) i^"""^ / Pm{R)dR . 

-oo J —oo ih m=l J ~oo 

(12) 

The tree geometry has the pathology that the number of sites near the boundary is 
proportional to the total number of sites in the tree. Surface effects are therefore strong. 
One can see this in (|T^) where Pi{R) is multiplied by the largest power of K. To eliminate 
these boundary effects and obtain an expression for the behavior of the bulk, one can 
use a subtraction procedure. Letting F„ be an extensive quantity in a tree of depth n, 
the quantity {Fn — KFn-ij/'i is the value of F per site far from the boundary when n 
is large [16]. To see this, note that the number of sites m steps from the boundary in a 
tree of depth n equals K times the number of sites m steps from the boundary in a tree 
of depth n — 1. The contributions to F from sites m steps from the boundary are thus 
cancelled in the subtraction. Under this subtraction, the number of sites we are left with 
is (K"+i + JsT" - 2)/{K - 1) - Ki^K"^ + iT""^ - 2)/{K - 1) = 2 sites which are far from 
the boundary. Thus dividing the difference F„ — KFn-i by two gives the value of F per 
site far from the boundary. Applying this subtraction to ([T2|) gives the average integrated 
density of states {uj{E)) per site far from the boundary in the limit n ^ oo 

{u{E)) = — / p,UR) dR / PrUR') dR' e(- - R') 

K+1 /-o 

+ / Preal(i?) dR . (13) 

Z J — oo 

For values of E and A such that there exists a complex fixed distribution Pcompiex(-R), one 
can use the relation (||) between Preai and Pcompicx to express {uj{E)) in terms of this complex 
distribution (see the appendix) 

{iuiE)) = -Imj^— J Pcomplex(i?) dR J Pcomplex(i?') dR' log(P' - -) 

- / Pcomplcx(i?) dR log(P)} (14) 
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The conducting properties: 

We turn now to the relation between P{R) and the conducting properties of the system. 
Imagine the situation shown in figure 2. Attach a wire to each boundary site of a branch 
of a tree as shown in figure 2. Suppose then that one sends a plane wave in at the left 
and allows it to scatter off the tree. For a branch of finite depth, some of the incoming 
wave will be reflected and some will propagate through the branch to the wires on the 
right. Of interest is what happens when the depth of the branch becomes large. There 
are two possibilities: either the wave is entirely reflected or some of the wave propagates 
through the tree into the wires on the right. In the former case, either there is a gap in 
the energy spectrum or the wavefunctions of the tree are localised, in the latter case the 
wavefunctions are extended. Now it is known that for a potential with unbounded support 
(e.g. a Gaussian) [18-20] there are states at all energies. In this reflection amplitude 

of modulus 1 (complete reflection of the wave) implies that the states at that energy are 
localised. Therefore, to determine whether at a given value of E and A the particle is 
localised, it suffices to compute the reflection amplitude in this experiment. 

To relate the reflection amplitude to the Ricatti variables, one needs to analyse what 
happens at the boundary. First consider the Schrodinger equation for a site a adjacent to 
the right boundary 



where b and c are the neighbors of a on the tree and in the wire respectively (flgure 2). In 
the wire the wavefunction is a plane wave ip = e^^^ going to the right (corresponding to no 
incoming flux from right inflnity) so that tpc/'ipa = The Ricatti variable Ra = ipb/'^Pa 
thus equals 



The wave vector k can be adjusted by putting a uniform potential on the wire. We note 
that the Ra for all sites adjacent to the boundary have a negative imaginary part and it is 



Starting with the boundary Ra given by ([T6|) , we now iterate the recursion (H) up to 
the flrst R in the wire on the left side of the tree (flgure 2). At this boundary we have 
both an incoming and an outgoing wave, so that the wave function on the left wire has the 
form ip = e**^^ + re~^^^ where r is the reflection amplitude. The reflection amplitude r is 
therefore related to the Ricatti variable R^ by 



'a 



(15) 



Ra = E-XVa- e 



(16) 



easy to check that the iteration (|^) preserves this property. 




—ik 



+ re 



ik 



1 + r 



(17) 
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or equivalently, 

Rd - e^'^ 

It is clear that if Rd is real, the numerator and the denominator are complex conjugates 
and so |r| = 1. It is also easy to see that if Rd has a negative imaginary part, |r| < 1. Thus 
if the imaginary parts of the initial Ri iterate to zero, the particle is totally reflected and if 
they do not, the particle has a non-zero transmission coefficient. To find the mobility edge 
it therefore suffices to study the complex fixed distribution Pcompiex(-R) and to determine 
at which values of E and A the probability of finding a non-zero imaginary part of R first 
vanishes. 



Remark: An idealized scattering situation like the one shown in figure 2 can be used 
in other cases including finite dimensional lattices to decide whether a system is conduct- 
ing or insulating. One could have for example an incoming wave on a wire attached to an 
internal site and outgoing waves at the boundary. A straightforward calculation but which 
would require new notations (that we will not present here) would show that a refiec tion 
amplitude smaller than 1 in the scattering situation is equivalent to the Green's function 
having a non- vanishing imaginary part when the energy E tends to the real axis [25,26 and 
references therein]. 

For arbitrary A, E and piV), one does not know how to find the fixed distributions 
that solve @) and in the following sections we will show how certain quantities can be 
expanded in the limit of weak disorder (A small). An exception is when the distribution of 
the potential is Cauchy [5] 

Using the fact that sums of Cauchy random variables are also Cauchy distributed and that 
the inverse of a Cauchy variable is Cauchy, it is easy to obtain the exact form of Preai(-R): 

^-i(i?) = - (20) 

7C [R — ay + 0^ 

where the parameters a and h are the values of the attractive fixed point of the following 
two dimensional map 

E~K-^ ■ h^^, = \ + K ^" 



a. 



Using (p!3D one then finds the following closed expression for the integrated density of states 
{u{E)) 

{uj{E)) = — tan-M- tan-M — . 21 
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Unfortunately, no one has to date been able to obtain an exact expression for the complex 
fixed distribution, even in the case of a Cauchy distributed potential. So conduction 
properties such as the location of the mobility edge are not known exactly even for Cauchy 
disorder. 



3. The pure system 

Before proceeding to the weak disorder expansions, let us consider the case of no disorder 
(A = 0). For zero disorder and for our choices of boundary conditions, all the Ri for the 
boundary sites are equal, and so the initial distribution Pi{R) is a 5 function. In the 
absence of disorder, all the Pk{R) computed from Pi{R) through the recursion (H) are also 
S functions concentrated at some value 

Pk{R) = S{R-Ak) (22) 

(23) 

E\ > 2\/K and complex 



where the A^ satisfy the following recursion 

Ak+i = E - — . 

This recursion for A^ has two fixed points which are real when 
conjugate when —2\fK < E < 2\/K. 



First, if 



E > 2\fK (24) 

the sequence A^ given by the recursion (^) always converges to the real fixed point A 
given by 

A = ^±^|IHM. (25) 

This implies that in the scattering situation described in section 2, the initial complex R^ 
become real under iteration, and hence the wave is completely reflected. Furthermore, the 
integrated density of states is zero for this range of energy. To see this, note that with 
the boundary condition (^ (i.e. Ri = E), appropriate for the calculation of the density of 
states, Ai > A2 > ■ ■ ■ > An ■ ■ ■ > A, and so all the R are positive. Moreover E — {K+1)/An 
is also positive. It then follows from expressions (0) or (|13D that the integrated density of 
states is zero for the range of energy (|2^), meaning that even for a finite tree, there is no 



eigenenergy greater than 2\'K and so the range of energy ( p4[ ) is outside the band. The 
case E < —2\/K is obviously symmetric. 
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On the other hand, if 

-2Vk<E<2Vk (26) 

there exists a complex fixed point of (^) 

, E - W4:K - , , 

A = . (27) 

meaning that there is a (5-distribution concentrated at this point that solves (0). 

If one starts with a real Ai as in (|^) to compute the density of states, the sequence Aj^ 
does not converge. Letting the energy E = 2\fK cos9, one finds 

/ — sinfn^ + 9) , ^ 

A^ = Vk ) ' . 28 
sm(n6') 

The integrated density of states for a tree of depth n then follows from (PI): 

^^(^) ^ e f """^ - """"("^ ^''A + iK + i)± e f- "-'A . (29) 

\ sin(n6' + ^) J \ sin(m^) / 

The spectrum consists of a finite number of eigenvalues (as it should for a finite system) 
with huge degeneracies which refiect the symmetries of the tree. 

If one starts with a complex Ai, the sequence (^) does not converge either. It is 
however easy to show that 



Ak-^e^' ^-2(.-i)..Az.^^. (30) 



Ak - \fKe-^^ - v^e-*« ' 

From this explicit expression, we see that if Ax is complex. A], remains complex, and so 
the system is conducting. 

Remark: tt is a property particular to the pure system that when the i?j at the boundary 
are all equal, the sequence An does not converge, and thus the sequence Pk(R) has no limit. 
As soon as one introduces disorder (A 7^ 0), there are limiting distributions Preai(-R) and 
-Pcompiex(-R), which Satisfy the fixed point (^. These distributions, in the limit A ^ 0, also 
satisfy the fixed point equation (^. It is possible to find these distributions. The complex 
distribution is a 5-function concentrated at A given by ( pTj) and the real distribution is 
given by the relation (^ between Preai(-R) and -Pcompiex(-R) 



- 2^ R^-ER + K ■ 



It is interesting to notice that this fixed distribution (|3lD is the invariant measure of the map 
(P^). Using this fixed distribution, one obtains the following expression for the integrated 
density of states uj{E) per site far from the boundary 

K-1, _^/K + l 



u{E) = -^e - tan-^ (^-^ tan ) . (32) 
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4. Weak disorder expansion inside the band 



We turn now to the weak disorder expansions. As the properties of the pure system are 
qualitatively different in the band and outside the band, the weak disorder expansion 
requires rather different techniques in these two cases. We will see that the neighborhood 
of the band edge must also be treated separately. 

In this section, we explain the small A expansion for energies E inside the band of the 
pure system < E < 2^/K). In this energy range, we saw that the complex fixed 

distribution solution of (|^) is a 5-function concentrated at the complex number A 

A^^-'^f^. (33) 

To obtain the weak disorder expansion we assume that for A small, the variables R which 
appear in the recursion (0) have small fluctuations around A (see figure 3 and remark 1 at 
the end of this section), so that the distribution Pcompicx(-R) remains concentrated around 
A. In accordance with this assumption [16], it is convenient to make the following change 
of variables ^ 

= 1 + XB, + X^Q + ... ■ ^^^^ 
Here Bi,Ci,-- - are fluctuating quantities, the distributions of which are assumed to 
independent of A. Inserting into (0) and equating terms order by order in A, one 
finds that the fluctuating parts Sj, Cj, ■ ■ ■ of i?j satisfy the following recursions 

AB, = V+jj:B, ■ AQ-AB^ = jj:C, (35) 

and so on. Using these equations, it is easy to compute the moments of the fluctuating 
variables Bi, Ci, ■ ■ ■ in terms of K, E and the moments (V^^) of the random potential V 
which we assume to be finite. For example 

Using this weak disorder expansion one can then calculate perturbatively any quantity that 
can be expressed as a function of the Rt. For example the real fixed point can be obtained 
from d^) 

The integrated density of states can also be calculated using ([1^) which relates u^E) to 
the complex fixed point. Substituting the perturbative expansion (Rl) of R around the 
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complex fixed point into (p!4[) , expanding in A and calculating the necessary moments we 
find to order 

ujiE) = — Im \K\ogA -^-^^^logiA' - 1) - X'{V')--^^ - A^(V^=^), 



TT [ ° 2 ' ' '2(^2-1)2 ^ '3(yl2- 

Higher order terms in A are straightforward to calculate but in practice the algebra required 
to work out the necessary moments and correlations becomes rather involved. 

Conduction properties can be calculated in a similar way. A question of basic interest 
is whether the wave function is localised or extended. From the discussion of section 2 
we know that the wavefunction is extended if there exists a complex fixed distribution 
concentrated in the lower half plane. One way of determining whether such a complex 
fixed distribution exists is to calculate the magnitude of the reflection amplitude and see 
whether it is less than (extended) or equal to (localised) 1. The existence of the complex 
fixed distribution can be determined more easily however by considering (Imi?) . Since the 
complex fixed distribution is concentrated in the lower half plane, a necessary and sufficient 
condition for it to exist is for (ImR) to be non-zero. Perturbatively we have 

(ImR) = (Im {A{1 - XB - X^C + X^B^..]) (39) 



Remark 1: Equation ( |39D indicates that for energies E inside the band of the pure system, 
weak disorder leaves the eigenfunctions extended. This result, obtained for arbitrary K 
including K = 1 where the tree becomes one dimensional, apparently contradicts the well 
known fact that weak disorder localises all the eigenfunctions in one dimension [27]. The 
reason the above calculation does not apply when i^' = 1 is that the assumption that the 
iteration @) converges for weak disorder to a fixed complex distribution close to the fixed 
point A is not valid. This can be easily seen in the weak disorder expansion. Assume that 
some Ri obtained after n iterations has the form Ri = A/{1 + Asf"^ + ■■■). Clearly one has 

5r^ = ^ + ]^Ei?r^ (40) 
This implies the following recursion for (|i?*^"^|2) 

(|5W|2) = il_i + iL^ \A (41) 

^ ' K K ^ ^ 
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where we have used the equahty \A\ = yK. We see from (^) that for K > 1 the fluctu- 
ations saturate as n increases, meaning that an initial distribution concentrated 
in the neighborhood of A remains concentrated near A. On the other hand, for K = 1, 
the fluctuations grow linearly with n so that the recursion does not converge to a fixed 
distribution close to A. 

Remark 2: The case K = 1 (one dimension) is special in other respects. In the band, 
at energy E = 2^/Kcos6, the fixed point A is A = V^e~*^. Hence for K = 1, there are 



a large number of values of 9 at which the weak disorder expansion (38) contains small 
denominators [17,28-31]. This is not the case for K > 1 where only the band edges (^ = 
and 6 = 7t) seem to give rise to small denominators. 

5. Near the Band Edge 

As E approaches the band edge of the pure system 2\/K from below at fixed A, terms 
like {A"^ — K)~^ diverge, and the above perturbation theory breaks down. Perturbative 
information about the complex fixed point can, however, still be obtained by choosing the 
appropriate scaling of the energy with A. Let 

i2 



E = VK{2 - aX^) (42) 

where a is of order one. For these energies, A is given by v^(l-iAv/a+0(A2)). Substituing 
this in the expansions (|38| ) and ( ^9]) one finds 

- ' 3^k^^ r ^ ^ 8(ir-l)U + ■ ■ -j ^''^ 



and 



(44) 



From this we see that the leading term in a small A expansion with a fixed as A — > 
corresponds to resumming an infinite series in the original expansions ( p8| ) and (|39D. 
For the range of energies (|42|), let us make the following change of variable 



R^ = ^ . (45) 

Substituing this expression in (|), we see that the 0i are of order A. Expanding both sides 
of (0) yields 

ct>, + ct>i + ctfi + ... = --^ - + E*^. (46) 
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The moments of can then be calculated by taking successive integer powers of this 
equation. To lowest order one obtains 

{<P') = -aX' ; ('^)' = -^'(«+i^) ' (# = 0(A^) (47) 

and so forth. Just as in the last section we can now calculate the average of functions with 
respect to the complex fixed point perturbatively. For example, one finds using (IH) for 
the integrated density of states 

and {uj{E)) = to order for a < -{V^)/{K - 1). Both (WTi) and agree with(p|) 



and (|39D when a becomes large as it should and (^71) andpg) are expected to represent the 
ressumed series of the most singular terms in the expansions ( p8D and (|39|). 

One might wonder how this result is modified by higher order terms. If one decides to 
calculate the integrated density of states or any other quantity to a given order p in A, one 
needs to calculate the first p moments to order A^. To do so, one takes the successive 
powers of (^) and one averages. This gives relationships between the first p moments of 0. 
If one uses these relations to express all the moments in terms of the first moment (0), one 
ends up with a polynomial equation in of degree p which has real coefficients depending 
on A, a and the moments of the potential V. As one varies a (i.e. the energy E), one finds 
that two complex roots of this polynomial become real at a certain critical value Ec which 
to fourth order in A reads: 



K-1 {K-iy (is: -1)3 A{K-1)^ 

If one pushes this expansion to higher order in A, one would find higher order corrections 
to Ec but for E > E^, (0) as well as all the higher moments of would be real, and this 



implies through (|T^) tha t id{E) = for E > Ec to any order in A. 

This result seems to contradict the well known fact that the density of states never 
vanishes for distributions of the random potential with unbounded support (such as the 
Gaussian distribution) [18-20] . Our result can be reconciled with a non-zero density of states 
at all energies E only if there are small non-perturbative contributions to the density of 
states above Ec which vanish to any order in A. We discuss this point in the next section. 

If perturbation theory could be trusted, it would follow from the discussion of section 
two that the particle is localised for energies above Ec since (0) is real. The perturbative 
calculation therefore predicts a mobility edge at Ec given by (^9]). This result could be 
modified by non-perturbative effects; unfortunately, unlike the case for the density of states, 
we have not been able to discover what these effects might be or how to calculate them. 
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6. Beyond the Band Edge 



In this section we consider energies larger than lyK and outside the scahng regime 
discussed in the last section. In this range of energy, the i?j converge to the real fixed 
point A given by (^) when A = 0. For small A, one expects the i?, to fiuctuate about 
A. Expression (|13D for the density of states requires however the knowledge of Preai(-R) for 
values of R far from A. This is a large deviation problem since untypical values of R are 
produced by untypical values of the potential. We are going to show that the real fixed 
distribution Preai(-R) can be calculated by a saddle point method. We have also tried to 
calculate the complex fixed distribution in this range of energies by a similar approach but 
failed and we are unable to tell whether the fixed complex distribution simply does not 
exist or whether we have just not been able to find it. 

Since the regions of R which contribute to the density of states corresponds to untypical 
values of i?, the shape of Preai(-R) in these regions depends strongly on the shape of p(V^). 
In what follows, we will take piy^ to be Gaussian 

P(^) = ^exp(^- (50) 

The equation to be solved for Pt-csi{R) is is then 

P{R) = / ^ exp (^] n dR^ P{R.) 6(R-E + XV + j:j-] (51) 



^ \ J i=l \ i=l Ri 

Now suppose that A is small and that R has fiuctuations of order A around A. If we write 
R = A + X(f), we get from (g) 

= -Ar + -1^0, . (52) 

Since V is Gaussian distributed, it follows that (p is also Gaussian with zero mean and 
(02^ = {V^)/{1 - K/A^). In terms of the variable P, 

valid for R — Aoi order A. Given this expression for P{R) valid for small fluctuations away 
from A, it is natural to look for the solution of (ISTF) of the form 



P{R) = Q{R) exp iJ^^ (54) 
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where F{R) is independent of A. Plugging this into (^T]) and integrating over V gives a 
saddle point equation for F{R) 

The saddle point is given by the Ri that solve 

m) = -;^(^-^ + E^) • (56) 



Finding the solution F(R) of (|55| - ^6D is non trivial. If however the saddle point is symmetric, 



i?i = i?2 = ■ ■ ■ = Rk, the function F{R) satisfies 



F{R) =meixl- 1(e - R- ^] + KF{Ri)\ . (57) 



[ 2 V RiJ 

Even under this assumption, we could not find a closed expression for F{R) in terms of 
elementary functions. One can however show that a solution of ( |57| ) is given by 

F(R) = lim GJR) (58) 

n— >oo 

where the Gn{R) are defined by 

G„(i?)= max \-ltK'(x,-E+^) ] (59) 

with Xq = R and Xn+i = A. 

It is easy to show that (|58| ) and (|59| ) solve (|57|). First, it is clear that Gn{R) < for 
all n and R. Second, it is easy to check that Gn{R) increases with n. This can be seen 
by choosing the same set Xi,X2, ■ ■ ■ Xn for Gn and G„+i with X„+i = X„+2 = A. Since 
the Gn increase and are all negative, the limit (|58D which defines F{R) exists. That F{R) 
satisfies { ^7\ ) then follows from the observation that 

G'„+i(i?)=niax|-^(^E-i?-^) +KGn{Ri)^ (60) 

We tried to check the assumption that the saddle point equation (^) has a symmetric 
saddle point by solving this equation numerically for K = 2 and E = 3. . The values of -Ri 
and i?2 which give the saddle point are shown in figure 4 and the numerical solution for 
F[R) is shown on figure 5. We see that the saddle point is symmetric for a range of values 
of R around A. Outside this range, however, the saddle point is no longer symmetric and 
F[R) seems to be constant. The rather complicated shape of F[R) makes the calculation 
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of the integrated density of states difficult, in particular because it is hard to tell if the two 
terms which appear in the expression (|13D of {uj{E)) are of the same order when A is small, 
and because if the two terms are of the same order, the prefactor Q{R) in (Q) would need 
to be calculated. In any case, one would find a non-zero integrated density of states {uj{E)) 
which would be exponentially small when A — 5> and which would vanish to al 1 orders in 
perturbation in A. 



7. Numerical attempt to determine the mobility edge 



The expression (^9]) obtained in section 5 predicts that the mobility edge tends to the 
band edge value E = 2^J~K as A 0. This disagrees with the result of Abou-Chacra and 
Thouless [5] that the mobility edge tends to = K + 1 as A gets small. 

That the perturbative expansion might be insufficient to predict the location of the 
mobility edge has been discussed above: the perturbative expansion predicts that for 
E > Ec the density of states is zero to all orders in A, whereas one knows that the density 
of states never vanishes for an unbounded distribution piV) of potentials. We tried to show 
in the previous section that a way of resolving this difficulty is to have non-perturbative 
contributions. The same could also happen for the reflection amplitude: one could have 
\r\ = 1 ioT E > Ec to all orders in A but with |r| 7^ 1 because of non-perturbative effects. 

One way of locating the mobility edge is to determine whether the distribution Pj.ea\{R) 
is stable against imaginary perturbations. We know that if the Ri are real at the boundary, 
they remain real under the iteration (H). Now let us add an infinitesimal imaginary 
component ie to the Ri at the boundary. After a finite number n of iterations of (^, the 
imaginary part of Ri is still infinitesimal and is proportional to e. Calling this imaginary 
part ieYi, we have 

y^=T. n ^ (61) 

paths J £ path j 

where the sum in (|6lD runs over all the K'^ paths of n steps from site i to the boundary 
of the tree. If one computes the Yi, either they go to zero as n increases, meaning that 
Preai(-R) is Stable or they grow with n, meaning that Preai(-R) is unstable. 



7. 1. Numerical approach 

We did not find an analytical way of determining whether the Yi grow or decay under 
the iteration procedure. So we had recourse to a MonteCarlo method. We represent the 
distribution Pi.eai(-R) by a sample of N points, where typically N = 100, 1000, 10000. This 
means that we have N values of the Ri and values of the Yi. At each elementary step. 
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we update one i chosen at random, by chosing K indices at random ji, ■ ■ -jx between 1 
and N and a random value of the potential Vi and we replace Ri and Yi by 



Ri 




(62) 



(63) 



We start with the of order 1 and we iterate this procedure until all the Y^ have become 
either very large 10^'' or very small 10"^''. The mobility edges estimated by this procedure 
are shown on figure 6 and 7 for the Gaussian and Cauchy distributions respectively. The 
points for different values of A are obtained for different samples, and so the roughness of 
the curves indicates the statistical errors. 

We expect this procedure to give the true mobility edge for infinitely large since 
for large the values of Rj^, ■ ■ ■ Rj^ are independent. We see that as increases the 
estimated mobility edge moves upward and to the right, and so it is not so easy to predict 
from these data accurate values of the large N limit. For a Gaussian potential and for 
small A, the results shown in figure 7b seem to agree well with the expression (|49|) . Notice 
however that the agreement becomes worse as increases. This means that even though 
the mobility edge in the limit A — > starts at 1^~K and not + 1 in our simulations, we 
cannot conclude from the data that the mobility edge really starts at 2\/K when A — > 
because the limits N ^ oo and A — > may not commute. 

7. 2. The case of independent R 

In order to test the validity of this numerical approach, it is useful to try it in an exactly 
soluble case. If we suppose the Ri to be independent random variables distributed according 
to a given probability distribution Pi.eai(-R), one knows from the theory of directed polymers 
[32,33] on a tree the exact expression of the large n limit of logl^/n where the Yi are defined 
by (|6lD. To test the MonteCarlo procedure descibed above, we computed the "mobility 
edge" by using (|63D with independent Ri chosen according to their exact probability 
distribution (^). Then using the known results from the problem of directed polymers 
one has for large n, 



The line in the plane E, A which separates the region of very large Yi and very small Yi can 
then be obtained exactly using the following expression for 




(64) 



{R-'') = - / 

TT J-oc 



dR 



{R-a)^ + h'^ cos(7r/3) 



(^)'-(2/^tan-0) (65) 
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The exact curve for independent Ri is shown in figures 6c and 7c together with the results 
of the MonteCarlo procedure. We see that the results seem to converge rather slowly as 
increases. Moreover, the limits N oo and A — * do not seem to commute, as the exact 
result tends to -ft^ + 1 as predicted by Abou-Chacra and Thouless [5] whereas the finite N 
results tend to the band edge value 

In summary we see that the MonteCarlo procedure described above can in principle be 
used to determine the mobility edge. However, = 10000 does not seem to be large enough 
to resolve numerically the question of where the mobility edge starts in the lambda 
limit. 

8. Conclusion 

We have shown that a great deal of information on the Anderson model on a tree is is 
contained in the random recursion (^. If we assume that under iteration of (^, the 
distribution of R converges to a fixed distribution P{R), the problem of knowing whether 
the wave functions are extended or localised reduces to the question of the existence of 
a complex fixed distribution Pcompiex(-R)- The equation (§) that P{R) satisfies is in fact 
equivalent to those already obtained by other methods [5-9]. We think however that our 
way of deriving these equations is more direct. 

The main result of this work was to show that one can expand quantites of interest 
like the density of states, the mobility edge or the refiection amplitude in powers of A. This 
approach is not however entirely satisfactory because we have not fully understood the non- 
perturbative effects. In particular, the existence and nature of non-perturbative corrections 
to the complex fixed distribution Pcompiex(-R) remains an open question, as, by consequence, 
does the position of the mobility edge in the limit A — 0. Both our perturbative expansion 
and our MonteCarlo simulations indicate that the mobility edge starts at the band edge of 
the pure system, but non-perturbative effects could change the former prediction, and as 
is the case for the independent R of section 7, the limit N ^ oo and the limit A — > may 
not commute in the latter. 

We think several important points deserve further consideration. 

First, it would be nice to be able to mathematically prove that the sequence of Pn{R) 
converges and to know under what conditions a complex fixed distribution exists. This 
problem is not easy because, as we discussed above, in the pure case, the sequence does 
not converge for energies inside the band and somehow it is the effect of a weak disorder 
which makes the distribution concentrate aroud the complex fixed point. 
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Second, it would be interesting to develop a non-perturbative approach, especially 
for the complex distribution, in order to compute at least for small A the shape of the 
mobility edge. Despite our efforts, we were unable to find a method allowing us to describe 
Pcompiex(-R) for Small A in the range of energies \E\ > 2y^. 

One could also try to use the recursion (H) to calculate other quantities which play an 
important role in the localisation problem, such as the inverse participation ratio. 

Last, we think that the MonteCarlo procedure described in section 7 could be used to 
accurately determine the position of the mobility edge by increasing to 1 or 10 million, 
though this would require a rather serious numerical effort. 



Appendix 



In this appendix we derive the relation (|^) between the real and complex fixed point 
distributions and the expression (|14D of the integrated density of states from the expression 
(|l^). To do so let us asume that we have a sequence of complex distributions -Pcompicxl'"^^'^) 
concentrated on the lower half plane which satisfy the recursion (|^)). To each of these 
distributions we associate a real distribution Q"'{R) by 

We are going to show that the sequence of Q"'{R) thus defined also satisfies (H). 

From the recursion relation for the complex distribution -Pcompiex' '^^^ rewrite 
lI]) as 

/oo rco ^ roc roc 1 S /" 
dr ds II dn / ds,-— — — / p{V)dV (A2) 
-oo JO ~S{J-oo Jo IT [K — r)^ + S'^ J 



^complex(r. - 5{r-E + XV + Y. 12^72 ) ^(^ - E 72^72 ) 

i=l ' i ~^ 1=1 ^ ' i 

It is useful now to note two properties of Cauchy distributed random variables. A Cauchy 
distribution is a probability distribution C(x; a, h) defined by 

1 h 

C{x] a, h) = — .2 I A2 (^3) 

Ti [x — aV + 
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where a is real and b is real and positive. Let xi and X2 be Cauchy distributed 
random variables with distributions C{xi] ai,bi) and C(x2;a2,&2) respectively. Then 
the sum x = Xi + X2 is Cauchy with distribution C{x;ai + 02, &i + 62)- Similarly, if 
X is distributed according to C(x;a,6), the inverse y = ^ is Cauchy with distribution 
C{y; a/ {a? + 6^), h/ (a^ + 6^)). Hence, if the real variables Rk are Cauchy distributed 

with distributions C{Ri] rj, Si), it follows for fixed V that the real variable R defined by 

k 1 

R = E-\V-Y.^ (A4) 

i=i 

is distributed according to C{R;r,s) where r = E — XV — YA=i^i/ {'^f + sj) and s = 
J2i=i Si/ (si + Written out as an equation, this reads 

/oo ^ 1 e . ^ 1 



If we substitute (|A5|) into ([A^), and use ( [All ) namely that 

-00 JO n [R — TiY + sf ^ 

it follows that the Qn{R) also obey (^. Assuming that in the limit n ^ 00 the Q"'{R) 
converge establishes (|). 



Let us now see how (0) can be obtained from (0). The expression (|T3|) can be rewritten 

as 

2 f _ 1 /-oo /-oo 2 

(^(^)) = ^_ / P,,,i(i?) di? / Prcai(i?') rfi?' log(/?' - - -) 

TT L / J-00 J-00 K — It 

^^^+1 ^^^^^^^^ log(i?-«e)j (AT) 
z J— 00 J 

where — is an infinitesimal imaginary part. One can then replace PreaK-R) in this 
expression by (|^). Using the residue theorem to do the integral over the real variable 
R gives (0). 
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Figure Captions 

Figure 1: A Cayley tree with K = 2 and deptli n=4. 

Figure 2: A branch of a Cayley tree with wires attached to the boundaries to test the 
conducting properties. In the scattering situation, a plane wave is sent in from the 
left and allowed to scatter off the tree. 

Figure 3: A set of typical R obtained by the Monte-Carlo version of the recursion (0) 
explained in section 7 for = 1000. Here K = 4, E = 2.1, and the complex fixed 
point solution of (pT]), represented by a diamond, lies a.t A = 1.05 — 1.702i. The 
ellipse represents the set of points to which the initial R = 2 — i would map in the 
absence of disorder. For small disorder, A = .1 (fig 3a), the R quickly move off the 
ellipse and concentrate themselves around the complex fixed point. As the disorder 
increases, A = 1. (fig 3b) and A = 5. (fig 3c), the R spread out. Finally, for very large 
disorder, A = 30. (fig 3d), the R quickly become real, and the particle is localised. 

Figure 4: The solutions Ri{R) and R2{R) of the saddle point equation (|55|) for = 2 
and E = 3. as a function of R. Near the stable fixed point A, the saddle point is 
symmetric: Ri{R) = R2{R), while for larger and smaller values of R the symmetry 
is broken, Ri{R) ^ R2{R). 

Figure 5: Numerical solution of the saddle point equation ( ^5]) for E{R) with K = 2 and 
E = 3. . The function F{R) seems to be constant in the range where the symmetry 
between Ri{R) and R2{R) is broken. 

Figure 6a: Monte Carlo determination as described in section 7 of the mobility edge 
for Cauchy distributed disorder. Here K = 4. The number A^ of different R is 
100, 1000, 10000. Figure 6b: Same as figure 6a for the Gaussian disorder. Figure 6c: 
Same as figure 6a for independent R distributed according to the Cauchy distribution 

m 

Figures 7a,7b,7c: Enlargements of figures 6a, 6b and 6c showing small values of A. Figure 
7b shows a rather nice agreement with the result (|49D . However as A^ increases the 
agreement seems to become worse indicating that the limits N —>■ oo and A — ^> do not 
commute. This is even more apparent for the independently distributed R where the 
exact result converges in the limit A^Oto-E' = A' + las predicted by Abou-Chacra 
and Thouless [5] whereas the finite A^ curves converge to the band edge 2y/K. 
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